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Abstract 



CD In the factorized picture of semi-inclusive hadronic processes the naive time reversal-odd parton distributions exist by virtue of the 



gauge link which renders it color gauge invariant. The link characterizes the dynamical effect of initial/final-state interactions of 



the active parton due soft gluon exchanges with the target remnant. Though these interactions are non-perturbative, studies of final- 
es state interaction have been approximated by perturbative one-gluon approximation in Abelian models. We include higher-order 

Q gluonic contributions from the gauge link by applying non-perturbative eikonal methods incorporating color degrees of freedom in 
a calculation of the Boer-Mulders function of the pion. Using this framework we explore under what conditions the Boer Mulders 
00 function can be described in terms of factorization of final state interactions and a spatial distribution in impact parameter space. 
^ \ 
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Introduction 



Over the past two decades the transverse partonic structure 
of hadrons has been the subject of a great deal of theoretical 
and experimental investigation. Central to these studies are 
j the early observations of large transverse single spin asymme- 
\q tries (TSSAs) in inclusive hadron production from protonpro- 
ton scattering over a wide range of beam energies QjJ, |2| |3l |4fl. 
i-H Recently TSSAs have been observed in lepton-hadron semi- 
inclusive deep inelastic scattering (SIDIS) at the COMPASS Ja] 
t-H |and HERMES 0] experiments and at Jefferson Lab fl, 
C7N as well as in inclusive production of pseu do-scalar mesons in 
O proton-proton collins at RHIC (HI [II [3 El- While the naive 
;> parton model predicts that transverse polarization effects are 



trivial in the helicity limit lfl4ll . it has been demonstrated 1 15 
[TH 12 , 3 that soft gluonic and fermionic pole contributions 
& to multiparton correlation functions result in non-trivial twist- 
three transverse polarization effects. In addition theoretical 
work on transversity [ 19, 20, 21] indicated that transverse po- 
larization effects can appear at leading twist. Two explanations 
to account for TSSAs in QCD have emerged which are based 
on the twist- three and twist- two 1E1II22II2I I2IIE5I1 ap- 

proaches. Recently, a coherent picture has emerged which de- 
scribes TSSAs in a kinematic regime where the two approaches 
are expected to have a common description B26ll27ll28ll29ll . 

In the factorized picture of SIDIS [24, 30] at small trans- 
verse momenta of the produced hadron Pj ~ kj « ^[Q^ 
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the Sivers effect describes a twist-two transverse target spin-5 j 
asymmetry through the "naive" time reversal odd (T-odd) struc- 
ture, Af(x,k T ) ~ S T (Px % T )fi T (x,i*) where Vfi 1 
is the hard scale, kj is the quark intrinsic transverse momen- 
tum, and P is the momentum of the target. For an unpolarized 
target with transversely polarized quarks-s^, the Boer-Mulders 
function is given by Ah(x, k T ) ~ s T ■ (P X kr)hf(x, k\). 
Dynamically, T-odd-PDFs emerge from the gauge link struc- 
ture of the multi-parton quark and/or gluon correlation func- 
tions 1 26[ 32,33, 34l which describe initial/final-state interac- 



tions (ISI/FSI) of the active parton via soft gluon exchanges 
with the target remnant. 

Many studies have been performed to model the T-odd PDFs 
in terms of the FSIs where soft gluon rescattering is approx- 
imated by pert urbative one-gluon exchange in Abelian mod- 
els JH HIS EE H [M 53, El |H B . We go beyond this 
approximation by applying non-perturbative eikonal methods 
to calculate higher-order gluonic contributions from the gauge 
link while also taking into account color. 

In the context of these higher order contributions we perform 
a quantitative study of approximate relations between TMDs 
and GPDs. In particular, we explore under what conditions the 
T-odd PDFs can be described via factorization of FSI and spa- 
tial distortion of impact parameter space PDFs l44ll . While such 
relations are fulfilled from lowest order contributions in field- 
theoretical spectator models I45[ 14611 a model-independent anal- 
ysis of generalized parton correlation functions (GPCFs) [47] 
indicates that the Sivers function and the helicity flip GPD E 
are projected from independent GPCFs. A similar result holds 
for the Boer-Mulders function for a spin zero target |5H]- How- 
ever for phenomenology, it is essentially unknown whether the 
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proposed factorization is a good approximation. Here we focus 
on the transverse structure of the pion in terms of the impact 
parameter GPD H", and the Boer Mulders function for which 
there are very few studies. Recent lattice calculations [49] indi- 
cate that the spatial asymmetry of transversely polarized quarks 
in the pion is quite similar in magnitude to that of quarks in 
the nucleon which lends supports the findings in [5jj|. Further 
understanding of the Boer-Mulders function for the pion may 
provide insight into the explanation of large cos 2<p ~ hf ® hf 
azimuthal asymmetry (AA) observed in unpolarized n~p Drell- 
Yan scattering Jm, 52, 53]. This work also has direct impact 
on studies of AAs and TSSAs in unpolarized and polarized ttN 
Drell-Yan experiments proposed by the COMPASS collabora- 
tion. In the latter case the TSSA is sensitive to the the nucleon's 
transversity through the convolution of hf ® hi . 



2. T-odd PDFs, Gluonic Poles and The Lensing Function 

The field-theoretical definition of transverse-momentum de- 
pendent (TMD) parton distributions in terms of hadronic matrix 
elements of quark operators serves as the starting point of our 
analysis. A classification of TMDs for a spin- 1/2 hadron with 
momentum P and spin S was presented in Refs. 024l l29l 15411 . 
In an analogous manner, it is straight-forward to obtain the 
TMDs for a spin-0 particle from the correlator for a pseudo- 
scalar target. One encounters two leading twist TMDs for a 
pion, the distribution for unpolarized quarks f\ and the distri- 
bution of transversely polarized quarks hf, the Boer-Mulders 
function. Adopting the infinite-momentum frame where the 
hadron moves relativistically along the positive z-axis such that 
the target momentum P has a large plus component P + and no 
transverse component we use the light cone components of a 



4-vector a* = 1 / V2(a° + a 3 ), 



(a ,a + ,a J_ ). The Boer- 



Mulders function, defined for SIDIS reads 



2e'Jkihi(x,jcl) = m n f — 

T T 1 T ' J 2{2n) 

x<P| qj(Q) [0 ; oon] ia i+ y 5 [^n + z T ; z] ?«(z) \P), (1) 



Z T ixP + ;r-ik T -z 
) 3 



where [x ; y] denotes a gauge link operator connecting the two 
locations x and y and the light-like vector = (1,0,0). Pos- 
sible complications with slightly off-light cone vectors as sug- 
gested in TMD factorization theorems rf3QL 55 1 are discussed 



below. Throughout this analysis we work in a covariant gauge 
where the transverse gauge link at light-cone infinity is negligi- 
ble. The gauge link in ([TJ is interpreted physically as FSIs of 
the active quark with the target remnants Q2l 13311 a nd is nec- 
essary for "'naive" time-reversal odd TMDs 1 22, 25, 31] to ex- 
ist 13311 . The Boer-Mulders function appears in the factorized 
description of semi-inclusive processes such as SIDIS or Drell- 
Yan SSSSSESEIEISIII in terms of the first 
£7-moment, 2m l Jif i \x) - j d 2 kxkjhf(x,kj). It was shown 
in Ref. lE6ll that the first &r-moment of the Boer-Mulders func- 
tion can be written in terms of a gluonic pole matrix element. 
Transforming the pion states in Eq. ([TJ into a mixed coordinate- 
momentum representation 11461 16211 results in an impact param- 



eter representation for the gluonic pole matrix element, 

(k T )(x) = m n hf\x) = f d 2 b T ^-e ixP+z ~ 

X<P + , 0V| q(Zl) hi ; 22] IXzi) cr i+ q(z 2 ) \P + , 0V>, (2) 

where the impact parameter bj is hidden in the arguments of 
the quark fields, z^ 2 = Tyn^ + bj and the 4-vector b^ T = 
(0, 0, b l T , b 2 T ). The operator /' originates from the time-reversal 
behavior of the ISIs/FSIs implemented by the gauge link op- 
erator in ([TJ) and is given in terms of the gluonic field strength 
tensor F 1 ™, 



2/ ! 'fe) 



dy [zi ; y]gF + '(y) [y ; z 2 ], 



(3) 



withy'" = y'n! 1 + b^. 

Turning our attention to GPDs of a pion, they are represented 
by an off-diagonal matrix element of a quark-quark operator 
defined on the light-cone |63, 64, 65], where "in"- and "out"- 
pion states are labeled by different incoming and outgoing pion 
momenta p and p' . One encounters two leading twist GPDs 
for a pion, a chirally-even GPD and the chiral odd GPDs 
H* B48I1 . We use the symmetric conventions for the kinematics 
for GPDs J6U], P = \(p + p') and A - p' - p. The skewness 
parameter £ is defined by A + = -2%P + , and t — A 2 . The impact 
parameter GPDs are obtained from the ordinary GPDs via a 
Fourier-transform of the transverse momentum transfer A7- at 
zero skewness £ = 0. The chirally-odd impact parameter GPD 
is expressed as 



dz 
2X2n)~ 



JxP + z 



(P + ,0 T \q(Zl)[Zi;Z2]a- +i q(Z2) \P + ,0 T ) 
2V T d 

= — ^-p-wrc*,^). 

m* db\ 



(4) 



( H" describes how transversely polarized quarks are distributed 
in a plane transverse to the direction of motion. This distri- 
bution function represents a transverse space distortion due to 
spin-orbit correlations l67ll . A comparison of the first 

moment of the Boer Mulders function © and the first deriva- 
tive of the impact parameter GPD < H", Eq. (@J, reveals that they 
differ by the operator I' which represents the FSIs. In various 
model calculations l62i ro8ll the FSIs are treated such 

that the two effects of a distortion of the transverse space par- 
ton distribution and the FSIs factorize resulting in the relation 



2m\h\ (X \x) 



I 



d 2 b T b T ■ J(x, b T ) —^(x, P T ), (5) 
dbi. 



where I is called the "quantum chromodynamic lensing func- 
tion" 1 62 ] . This factorization ([5]) doesn't hold in general [48, 
l69tl . On the other hand it is unknown how well Eq. ([5]) works as 
a quantitative and possibly phenomenological approximation. 
A phenomenological test of Eq. (O requires information on the 
parton distributions and H* (in principle measurable) and 
quantitative knowledge of the lensing function. In the follow- 
ing sections we estimate the size of the lensing function using 
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Figure 1: The amplitude W including FSIs between re-scattered eikonalized quark 
and antiquark. The FSIs are described by a non-perturbative scattering amplitude M that 
is calculated in a generalized ladder approximation. Gluon interactions as shown in the 
second diagram are not taken into account (see text). 



non-perturbative eikonal methods |70l 21] to calculate higher- 
order soft gluon contributions from the gauge link and study 
how these soft gluons impact Eq. (0. Up till now the relation 
(|5]l was used to predict the sign of T-odd TMDs in conjunc- 
tion with numbers for the u- and cZ-quark contributions to the 
anomalous magnetic moment of the nucleon and the assump- 
tion that FSIs are attractive 16^1 . We will also investigate the 
latter assumption. 

3. TMD-GPD Relation for a Pion 

We focus our attention on a pion in a valence quark configu- 
ration that one expects for relatively large Bjorken x. Working 
in the spectator framework Q35i |3a HI, 72, 7311 and inserting 
a complete set of states, 1 = £ t in the quark corre- 

lation function Eq. (Q]), we truncate this sum to an antiquark 
and neglect multi-particle intermediate states. The usefulness 
of this approach is twofold: First, we are able to improve on 
the one gluon exchange approximation for FSIs to studying T- 
odd PDFs by including higher order gluonic contributions and 
color degrees of freedom. Second we are able to explore to 
what extent transverse polarization effects due to T-odd PDFs 
can be described in terms of factorization of FSIs and a spatial 
distortion of impact parameter space including higher gluonic 
corrections [46, 48] with color. Thus, we express the pion Boer- 
Mulders function (0 in the following way 



8(2tt) 3 (1 -x)P 4 



Y Wicr i+ y 5 W, 



(6) 



with the matrix element W given by 

Wf(P, k; cr) = (P - k, cr, S\ [oon ; 0] afi gf(0) |P>. (7) 

where <x and 6 represent the helicity and color of the interme- 
diate spectator antiquark. We model (0 by the diagram shown 
in Fig. where the FSIs - generated by the gauge link in (0 
- are described by a non-perturbative amputated scattering am- 
plitude (M)^ with /3, a (y, 5) color indices of incoming and 



outgoing quark (antiquark). In the next section we calculate 
the scattering amplitude using non-perturbative eikonal meth- 
ods thereby considering a subclass of possible diagrams with 
interactions between quark and antiquark. We neglect classes 
of gluon exchanges in the second diagram in Fig.0represented 
by the red rungs since they would be attributed to the "inter- 
action" between the quark fields and the operator / in (0 and 
lead to terms which break the relation ©. We also neglect real 
gluon emission and (self)-interactions of quark and antiquark 
lines the second diagram in Fig. 0since they represent radiative 
corrections of the GPD and are effectively modeled in terms of 
spectator masses and a phenomenological vertex function. 

The pion-quark vertex is modeled with the interaction La- 
grangian 



L 



(8) 



where we allow the coupling g„ to depend on the momentum of 
the active quark in order to take into account the compositeness 
of the hadron and to suppress large quark virtualities 11421 143 , 
7311 . Applying the Feynman rules we obtain an expression for 
the matrix element W in (0 from the first diagram in Fig.Q] 



W%(P,k)= ~ IT 



S^gAk 2 ) 



[(t+m q )v(P s ,cry}. 
k 2 —m 2 +i0 



r dS 

J (2« 



g„ ((P-q) 2 ) [(r-4+m q )y 5 (4-m s ) (M)$ (q, P s )v(P s , cr)]. 



[n ■ (P-q) + iO] [(P-q) 2 -m 2 +iO\ [q 2 -m 2 +i0] 



(9) 



where P s = P - k is the spectator momentum. The first term 
in (0 represents the contribution without FSIs while the sec- 
ond term corresponds to the first diagram in Fig. We then 
express the FSIs through the amputated quark - antiquark scat- 
tering amplitude M. Here both incoming quark and antiquark 
are subject to the eikonal approximation (see, e.g. B74I1 and ref- 
erences therein). While the active quark undergoes a natural 
eikonalization for a massless fermion since it represents the 
gauge link contribution, the eikonalization for a massive spec- 
tator fermion is a simplification that is justified by the phys- 
ical picture of partons in an infinite momentum frame. The 
eikonalization of a massive fermion can be traced back to the 
Nordsieck-Bloch approximation lp75ll which describes a highly 
energetic helicity conserving fermion undergoing multiple scat- 
tering with very small momentum transfer. In this approxi- 
mation the Dirac vertex structure u(p\)y^u(p2) ~ pp/m = 
where (p\ + pi)j2 = p. For a massive anti-fermion one iden- 
tifies the velocity v 11 - —p 11 / m, and the numerator of a fermion 
propagator becomes i(-p + m) — > i(—v ■ p + m). 

We proceed by performing a contour-integration of the light- 
cone loop-momentum q~ in Eq. (0 where we consider poles 
which originate from the denominators in (0. This assumes 
that the scattering amplitude M does not contain poles in q~ 
and the integrand is well behaved on the contour in q~. Be- 
fore we proceed, it is important to point out that one-loop cal- 
culations of T-odd functions were performed in a scalar di- 
quark model J33 ,35,38] and a quark target model IT76Tl where 
there are no contributions from a pole in q~ in the exchanged 



gluon propagator. This is one reason why a factorization of 
the form (0 is exact in the one gluon exchange approximation. 
One does not expect this feature to hold in higher order cal- 
culations. In fact when including axial-vector diquarks i42ll in 
the one-gluon exchange approximation a q~ pole contribution 
from the exchanged gluon exists which leads to light-cone di- 
vergences when q~ — > 00 and q + — > 0. An introduction of a 
slightly off-light-like vector v regulates this divergence result- 
ing in logarithmic dependence of the form, log(u + /iO [77]. 
Such logarithmic terms prevent a factorization of the form (0. 
Alternatively, one may introduce certain vertex function g„(k 2 ) 
by hand that suppress contributions from q~ - poles J42I1 . Per- 
forming the contour integration on q~ under these assumptions 
fixes the momentum q~ of the antiquark in the loop in (0 to 
q- = (q\ + m 2 s )l2q + . 

The eikonal propagator can be split into a real and imaginary 
part using l/(x+/0) = T(\ / 'x)—in6(x). It has been argued in [46] 
that only the imaginary part contributes to the relation (0 as 
it forces the antiquark momentum q to be on the mass shell. 
Thus, the imaginary part of the eikonal propagator corresponds 
to a cut of the first diagram in Fig. [TJ From the point of view 
of FSIs, the kinematical point q + — (1 - x)P + is the 'natural' 
choice for the plus component of the spectator. In the picture 
where one imagines the scattered quark and antiquark to move 
quasi-collinearly with respect to the target pion - backwards 
and forwards respectively - the quark and antiquark exchange 
soft gluons. Under these kinematic condition one would expect 
the FSIs to be dominated by the "small" transverse momenta 
of quark and antiquark rather than the "large" plus momenta. 
An integration over q + in (TTOb where contributions other than 
the pole term contribute include configurations where large mo- 
mentum is also transferred from quark to antiquark in the plus 
direction. Nevertheless the principle value does contribute to 
the integral ( TTOb which allows for such momentum configura- 
tions. While this effect is beyond the picture of FSIs from soft 
gluon exchange, we will consider this in a future publication. 
Proceeding with the picture of soft gluon exchange there is a 
clean separation of FSIs and spatial distortion of the parton dis- 
tribution in the transverse plane in the sense of 01. Using only 
the imaginary part of the eikonal propagator Eq. (0 reduces to 



[(jt + m q )y 5 v(P s ,cr)l 



f 



kj + m 2 

ad 



d 2 ftr {(P-q) 2 ) [(f-4 + m q ) 75 v(P s , a-)]. (MY (q; P s ) 



(2ff)2 



|g 2 + m 2 ] 



(10) 



We have introduced the notation M = m.,M/(2(l - x)P + ). 

Now we use ( fTOb to calculate the pion Boer-Mulders function 
via (0. Specifying the pion-quark-antiquark vertex function 



- gn (n _ 1)l A2 k 2_ A2+iQ . 



(ID 



where / is a homogeneous function of the quark virtuality, 
we choose it to be a Gaussian exp[-/l 2 |£ 2 |] in accordance with 



Ref. [42]. Inserting ( TTOb into (0 and a bit of algebra yields the 
following expression for the Boer-Mulders function 

<i>y T h\(x, P T ) = ^^L(xm s + (1 - x)m q ) ((1 - x)A 2 ) 2 '"' 



r d 2 q T 

J w 



d 2 p T 



e J j(q J T - p J T ) 



-j^ x (xm 2 s -xa-x)ml) --^(q^+p- T ) 



A 2 (x)]" [f T 



t) 2 (In) 2 TK ,r 
: (3[M eik ])^ (jcj + 4 T )((2n) 2 5 a Ptf 2) (p T + k T ) 



•A 2 (x)]" 



(nM eik ]fjjk T +p T )j, (12) 



with A (x) = - x(l - x)m 2 + (1 - x)A . Anticipating 
an eikonal form for the scattering amplitude M(x, kj,qj) — > 
M elk (|^7- + Afrl) that will be discussed in the next section we 
exploit this property to simplify the expression and show a re- 
lation to the chirally-odd GPD H\. Since GPDs are defined 
from collinear light-cone correlations functions gauge link con- 
tributions to GPDs don't lead to an observable effect. In fact, 
in light-cone gauge the corresponding contributions from the 
gauge link are re-shuffled into the gluon propagators [35] and 
they appear as gluon dressings of the tree-level contribution to 
GPDs. Thus one can consistently describe GPDs from tree- 
level diagrams in the spectator model where the effects of gluon 
dressings are effectively hidden in the model parameters. A cal- 
culation for the GPD H* for an antiquark spectator can be found 
in Q48D . It is easy to generalize it with a phenomenological ver- 
tex function (fTTT >. We obtain the following representation 

2 X 2«-' 



HUx,0, -Ai) 



2(2tt) 3 A 2 



(xm s + (1-x) nia) 



(l-x)A 2 



D 2 T +A 2 (x)) 



dip 
Jo 



dz 



z 2«-2g2i 2 A 2 e - 



1 -4z(l -z) 



Di 



Dl+A 2 (x) 



COS 2 If 



(13) 



where D T = i(l - x) Aj.. Performing a translation of the in- 
tegration variables in (fT2l i according to qj — » qj + kj and 
Pt pj+kj, a rotation of the form q' T —qj-pT,p' T = qr+Pr, 
weighting with a transverse quark vector k' T and integrating 
both sides over kj we find the relation 

\2\ 



m 2 X m (x) 



-I 



-^-^r q T ■ I(x, q T )H K , 



x, 0,- 



1 -x 



(14) 



The function /' can be expressed in terms of the real and imag- 
inary part of the scattering amplitude M, 



ft*. AO - ^ J^ 2 (2 P T-q T y (3[M^; ( |^|) 
((2tt) W 2) (/?r - qr) + (^[M eik ]^(|/ r - 



(15) 



In order to derive the relation (0 one transforms Eq. ( TT4T > into 
the impact parameter space via a Fourier transforms of the fol- 
lowing form, 



<W (x, b 



d 2 A T 



-iAj-Ej 



H^(x,0, -A 2 ). 



(16) 



The lensing function in the impact parameter space then reads, 



f(x,b T ) = i(l-x) 



J (2n 



\e l I'(x,q T ). 

Y 



(17) 



In the following section we will use a quark-antiquark scatter- 
ing amplitude computed in relativistic eikonal models as input 
for the lensing function < fT3T >. 

4. The Lensing Function in Relativistic Eikonal Model 

In order to calculate the 2 — > 2 scattering amplitude M 
(needed for ( flBl ) we use functional methods to incorporate the 
color degrees of freedom in the eikonal limit when soft gauge 
bosons couple to highly energetic matter particles on the light 
cone. It is non-trivial to extend the functional methods estab- 
lished in an Abelian to non-Abelian gauge theory such as QCD. 
Attempts in this direction were made in Refs. 11711 17811 . and only 
recently a fully Lorentz and gauge invariant treatment was pre- 
sented in Ref. [729. Here we outline the details of the functional 
approach as it pertains to implementing color structure to the 
scattering amplitude M and thereby the lensing function. We 
leave the details to a forthcoming publication [80]. 

Starting from the generating functional Z for QCD in a co- 
variant gauge, a quark antiquark 4-point function T can then be 
defined by functional derivatives with respect to quark sources 
which yields, 



T 2 - 



j DA e~i I^+W-d-Af) gTr In G" 1 [A]+Trln H~ l [A] G[A] Q [A] 

(18) 

The first exponential describes the gluonic part of the the- 
ory including self-interactions and the second exponential de- 
scribes internal closed quark and ghost loops. G, G are the 
non-perturbative quark- and antiquark-propagator determining 
the external legs of the 4-point function T, and H is the ghost 
propagator 18011. One imposes eikonal approximations on these 
propagators |p70l 17111 that simplify the computation of the path- 
integral. In an Abelian theory the eikonal approximation as dis- 
cussed in the previous section leads to a well-known eikonal 
representation [70], which was argued in B7ll 17811 to generalize 
to QCD in the following way, e.g. for a massless fermion 



Gf p (x,y\A) 



Jo 



t \ U-ig £ dpn-A a (y+pri)fV 
V /a, 



(19) 



where color is implemented by a path-ordered exponential in- 
dicated by the brackets (...) + and the color matrix f in the ex- 
ponential. 

Inserting the eikonal representation for the quark- and anti- 
quark propagator into Eq. (fT8l and implementing the general- 
ized ladder approximation one finds the color gauge invariant 
result corresponding to the picture of FSIs discussed in the pre- 
vious section, 

saS . -> (l-x)P + 
) (X. 13V + brft = ' 

/Sfi 



(M e %{x,Wr 



I 



d 2 z T e 



(20) 



f d N '^a f — e r ia - u ( e >x<l?T\)t-a\ ( it-u\ 
J J (2jt) n ?- 1 1 ,asK V 



In this expression, the - 1) dimensional integrals result 
from auxiliary fields a a (s) and u a (s) that were introduced in 
the functional formalism (see Ref. 17111 ) to separate the phys- 
ical gluon fields from the color matrices. The eikonal phase 
^(IfVl) in Eq. (1201 represents the arbitrary amount of soft gluon 
exchanges that are summed up into an exponential form and is 
expressed in terms of the gluon propagator in a covariant gauge, 



da I dpn^D^iz + an-pn), (21) 

oo U -co 



xQ 



where D denotes the gluon propagator, and g is the strong cou- 
pling. In this form the 4-velocity vector is expressed in terms 

(1 — x)P + — 

of the complementary light cone vector n where v = — 1 — - — n, 

with n-h-\ and n 2 = 0. One may choose n = (0, 1, Or). 
In Eq. d20T i we evaluate the color integral, 

f M af^af |^e--»(e^H^ L(e ,»^_^ 

(22) 

by deriving a power series representation for arbitrary A^ c . We 
expand the exponential exp[/yf • a] and rewrite the resulting 
factors as derivatives with respect to u. Then we perform in- 
tegrations by parts which reduces the a integral to a simple 5- 
function. This simplifies the M-integral where u is set to zero 
after differentiation We obtain 



d"(e' 7 > | 
du a K..du a -\u=o 



(23) 



Now we expand the remaining exponential in Eq. ( 1231 and note 
that one can write the set of partial derivatives with respect to 
u a ' as a sum over all permutations P n of the set {1, «}, which 
results in the power series representation for /, 



u<x) - hj^i h 



Ni-\ 

■Z2>" 

a„ = l P n 



<7] ^a /1 ^a Pn(i ) f<lp n (n) 



...r 



>V (24) 



This color factor matrix nicely illustrates the generalized ladder 
approximation. If only direct ladder gluons were considered 
the sum over permutations would become trivial in Eq. d24l and 
only terms (f 1 ...f"f" ...f') af} = C n F 5 af i with C F = would 
contribute. This constitutes the leading order in a large-A^ c ex- 
pansion while non-planar diagrams, i.e. crossed gluon graphs, 
are suppressed. For the leading contribution one may simply 
replace a — > Cf a s and work in an Abelian theory. In particu- 
lar, this replacement was suggested in perturbative model cal- 
culations EEO]. 

Since we take into account crossed gluons 
we have to sum over all permutations in (F24l . and such a re- 
placement is not possible. In an Abelian theory, the generating 
matrices t reduce to unity, t = 1, and since we have n! permu- 
tations of the set { 1 , n], we recover the well-known result for 
the Coulomb phase, 



*—i n\ 



(25) 



For the non-Abelian N c = 2 theory the generators are given by 
the Pauli matrices <x a = 2f. Instead of using the power series 
representation we can calculate the integral (1221 analytically by 
means of the relation (e'" ?) ff/J = ^cos(lfl) + ^dsin(^). 
We obtain a slightly different result compared to Ref. IT7lll for 
SU(2), 



*sugq,Xs s ( X X ■ X , , An. ■ X ,X X\\ 
Lp '(4) = So,, (cos - - - sin --1 + 1 (2 sin - + - cos -]] . 



As a check on our numerical and analytical approaches we nu- 
merically calculate the lowest coefficients in the power series 
d24l) . and they agree with the coefficients in an expansion in 
X of the analytical result ( |26*1 ). The disadvantage of using the 
power series representation ( l24t is apparent for numerical cal- 
culations since the number of operations grows with n\. That 
said, for SU(2) we calculated the first eight coefficients. For 
QCD, N c = 3, the generators t are given by the Gell-Mann ma- 
trices A a = 2f. Due to the difficulty of integrating over the Haar 
measure in Eq. d22l we put off the analytical treatment [80]. Us- 
ing the power series ( f24b we derive the following approximative 
color function for a — ^/4 

^tfap 1 = dafii-cia 2 + c 4 a 4 - c 6 a 6 - c 8 a 8 + ...), 

3[/^ (3) ](a) = M c i« - c 3fl 3 + c 5 a 5 - cW + ...), (27) 

with the numerical values c\ - 5.333, ci = 6.222, C3 = 3.951, 
c 4 = 1.934, c 5 = 0.680, c 6 = 0.198, c 7 = 0.047, c 8 = 0.00967. 
Working in coordinate space we express the lensing function 
directly in terms of the eikonal phase x defined in Eq. (fJTJ. 
Defining the eikonal amplitude as in section [3] where the real 
and imaginary part are 



3[M«j](£> = if 



d 2 Z^%[fafi(x(\ZT\))l 



(28) 



(29) 



we insert (f28]l and (1291 into the lensing function ( [TBI then trans- 
form it via ( TTTb into the impact parameter space. This yields a 
lensing function of the form, 



T'(x,b T ) 



C[|] 



d-x) Vrx'^ 



2N C 



\bi 



4 



(Tr3[/])' (|) + ^Tr [(3[/])' (|) (%[/]) (|) 



-I T r 

2 



(30) 



where ^' denotes the first derivative with respect to \zt\, and 
(3[/])' and (^R[/])' are the first derivatives of the real and 
imaginary parts of the color function /. Also, the eikonal phase 
is understood to be a function of |?rl/(l - x). Inserting ( l25l l 
into < f30l > results into the following expression for the lensing 
function in an Abelian U(l)-theory 



P um (xJ T ) = (l-x)^-x , ( 
4\b T \ 



\b 7 



1-x 



1 + cos^-( 



1-x 



(31) 



Similarly from (l26b we calculate the lensing function in an 
SU(2)-theory 



r su(2) (x, b T ) 



(l-x)b l T \b T \ 

— 5 — x(-, — ) 

l6\b T \ 1 " x 

>2 



(32) 



3(1 +cos^) + (^) -sin^ -sin^) 
v 4 U) 4\4 4/ 



(26) ^ or tne SU(3)-QCD case we use Eq. d27b . In Fig.|2]the function 



C[4] is plotted versus \. While the convergence of the power 
series is slightly better for SU(2) where the numerical result, 
calculated to eighth order, agrees with the analytical result up 
to 4 ~ 2, we can trust the numerical result computed with eight 
coefficients up to j ~ 1.5 for SU(3). 

At this point we discuss the eikonal phase x as defined in 
( |2"TT i which is determined by two quantities, the strong coupling 
g and the gluon propagator T>. One can write a general form for 
the gluon propagator in momentum space 



(33) 



where the gauge dependent part is in t>2- However, the gauge 
dependent part does not appear in the eikonal phase when in- 
serting Eq. ( 1331 ) into Eq. ( f2Tb because the eikonal vectors n and 
v ~ _ Lzgf ^ are light-like. Performing the integral yields the 
following expression for the eikonal phase 



(34) 



X(\zt\) = ^- dk T k T MtfMDii-ldt), 
2tt Jo 



where Jo is a Bessel function of the first kind. The gluon prop- 
agator represents all exponentiated gluons exchanged between 
the two eikonal lines in the generalized ladder approximation in 
Fig-Q] The couplings represent the strength of the quark (anti- 
quark) - gluon interaction in Fig.Q] 

As a check of the calculation we investigated the perturba- 
tive limit of our calculation. Assuming that the quark - gluon 
interaction g 2 is small and using perturbative gluon propagator 
in Feynman gauge for t>\ one can expand our non-perturbative 
result in Eq. (l30l to g 2 . The leading order corresponds to the 
result of the one-loop calculation of the Boer-Mulders function 
of Ref. ll48ll after additional eikonalization of the antiquark. 



5. Non-perturbative Quantities from the Dyson-Schwinger 
approach 

In order to obtain a numerical estimate for the eikonal phase, 
it is important to have a realistic estimate of the size of the QCD 
coupling g or a s = Since all the gluons exchanges between 
the eikonal lines are soft, the interactions take place at a soft 
scale. Thus we need to know the running of the strong coupling 
in the infrared limit. Inserting a perturbative gluon propaga- 
tor might not describe the gluon exchange realistically. One 
would expect that a non-perturbative gluon propagator would 




Fi glire 2 1 Left: C[j] of Eq. BOl as a function of the eikonal amplitude 4 , We compare the numerical result computed by means of Eq. 1241 up to the order n = 8 with the analytical 
result in Eq. 1331 for the S U(2) color case. The numerical and analytical result agree up to j ~ 2. For S{/(3), we compare the numerical results for the orders n = 7, 8. The results 
are reliable for | ~ 1.5. Center: The eikonal phase x° S (l-H) vs - \%t\ with input from Dyson-Schwinger equations at scales Aqcd = OGeV, 0.2GeV, 0.5GeV, 0.7 GeV. Right: The 
lensing function I'(x, bj) from Eq. 1301 for U(l), SU(2) and SU(3) for x = 0.2 at a scale Aqcd = 0.2 GeV. For comparison we also plot the perturbative result of Ref. 1 4-8] including the 
eikonalized antiquark spectator with an arbitrary value for the coupling, a = 0.3. 



be a better choice. The infrared behavior of both quantities, the 
running of the strong coupling and the non-perturbative gluon 
propagator, have been studied in the framework of the Dyson- 
Schwinger equations l82l 83, 84, l85ll and also in lattice(see 



e.g. isdn . One learns from such studies that the strong coupling 
has a value of about a s (0) ^ 2.972 in the infrared limit. In par- 
ticular in Ref. 1 82] fits were presented for the running coupling. 
Since we are merely interested in a numerical estimate of the 
lensing function we will apply the simplest form of the running 
coupling presented in [82], 



a s (0) 



In [e + a^/A 2 )" 2 + b^/A 2 )^] ' 



(35) 



The values for the fit parameters are A = 0.7 1 GeV, a i = 1 . 1 06, 
a 2 = 2.324, b x = 0.004 and b 2 = 3.169. These calculations 
were performed in Euclidean space where Landau gauge was 
applied, and agree reasonably well with each other. Because the 
light cone components in Eq. ( l34b are already integrated out and 
the remaining integration range is over a 2-dimensional trans- 
verse Euclidean space, and because the gauge dependent part of 
the gluon propagator does not contribute, it is natural to apply 
the Euclidean results in Landau gauge of the Dyson-Schwinger 
framework. One unique feature of Dyson-Schwinger studies of 
the gluon propagator is that it rises like (k 2 ) 2 "^ 1 in the infra-red 
limit with a universal coefficient k ^ 0.595. This makes it in- 
frared finite in contrast to the perturbative propagator. A fit to 
the results for the non-perturbative gluon propagator has been 
given in Ref. B0II1], 



Z(p 2 ,H 2 ) 



a Am 2 ) 



I I 2 \* I 2\ 2k 



(36) 



with the parameters c = 1.269, d = 2.105, and 5 = — Jj. 
These fits for the running coupling and the gluon propaga- 
tor merge with the spirit of the eikonal methods described 
above since closed fermion loops (quenched approximation) 
were neglected. By using the non-perturbative propagator ( l36l l. 



we partly reintroduce gluon self-interactions that were origi- 
nally neglected in the generalized ladder approximation. Ac- 
cording to Ref. [88]] the fitting functions Eqs. ( 133b and ( 1361 ) 
were adjusted to Dyson-Schwinger results obtained at a very 
large renormalization scale, the mass of the top quark, fi 2 = 
170 GeV 2 , which defines the normalization in d36b . Since the 
lensing function deals with soft physics, intuitively we pre- 
fer a much lower hadronic scale which sets the normalization, 
yU = Aqcd ~ 0.2 GeV. In the spirit of Sudakov form factors we 
also assume that the scale at which the gluons are exchanged 
is given by the transverse gluon momentum that we integrate 
over. In this way the running coupling serves as a vertex form 
factor that additional cuts off large gluon transverse momenta. 

Our ansatz for the eikonal phase given by Dyson-Schwinger 
quantities then reads, 



X DS (\zt\) 



= 2 dk T k T a s (k 2 T )J (\MkT)Z(k 2 T ,A 2 )lk\. (37) 
Jo 



The numerical result for this ansatz is shown in the center panel 
of Fig. [2] We plot this function for various scale Aqcd = OGeV, 
0.2 GeV, 0.5 GeV, 0.7 GeV. Although the choice of this scale 
is rather arbitrary we observe only a very mild dependence on 
this scale as long as it remains soft. We further observe that 
the phase doesn't exceed a value of 4 - 4.5 — > Xtmxl^ ~ 1.15. 
Thus this feature makes the application of the power series of 
the color function in SU(3) reliable since %/4 never exceeds 1 .5 
in the lensing function, Eq. ( f30b and in turn in the calculation 
of the Boer-Mulders function in Eq. (0. 

Finally, we insert our ansatz for the eikonal phase into the 
lensing functions (l30l > for a U(l), SU(2) and SU(3) color func- 
tion. We plot the results in Fig.[3]for a color function for U(l), 
SU(2), SU(3). While we observe that all lensing functions fall 
off at large transverse distances, they are quite different in size 
at small distances. However for each case, the all order calcula- 
tion sums up to an exponential of the eikonal phase where one 
observes oscillations from the Bessel function Jo of the first 
kind. Despite these oscillations the lensing function remains 
negative. 
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Figure 3: The pion Boer-Mulders function, xm~h l (a) vs. .v calculated by means of 
the relation to the chirally-odd GPD H" for a SU(3), SU(2), U(l) gauge theory. 



6. The Pion Boer-Mulders Function 

In this section we use the eikonal model for the lensing 
function together with the spectator model for the GPD H* to 
present predictions of the relation (0 for the first moment of the 
pion Boer-Mulders function 

We start by fixing the model parameters in ( TT3l . We en- 
counter six free model parameters m s , m q , A, A, g„ and n that we 
need to determine by fitting to pion data. In order to do so we 



determine the chiral-even GPD F* (for definition and notation 



see Ref. 114 810 in the spectator model, 



F\{x,0, -A 2 -) 



r*2n r*\ 

Jo Jo 



Sn _ e 2^ 2 A 2 



2(2n) 2 



2 ^ 



(1 - x)A 



ln-2 



(38) 



dz- 



z 2„-3 


, ry x(\-x)m^-2m s m q 


2A 2 2 ,+A 2 (x)) 
e (Hi 




I 1 4z(1 z) z5 2 + i(.v) COs2 ^ 


n 



When integrated over x, the GPD reduces to the pion form 
factor F" (Q 2 ) = —F" (Q 2 ). An experimental fit of the pion 
form factor to data is presented in Refs. lf89L l90ll . and up 
to Q 2 = 2.45 GeV 2 is displayed by the monopole formula 
FfniQ 2 ) = (1 + 1-85 Q 2 ) . This procedure is expected to pre- 
dict the f-dependence of the chirally-odd GPD H* reasonably 
well up to Q 2 = 2.45 GeV 2 . In order to fix the x-dependence 
of FF[ we fit the collinear limit F n x (x, 0, 0) to the valence quark 
distribution in a pion. A parameterization for F*(x) was given 
by GRV in Ref. 119 ill at a scale p 2 = 2 GeV 2 . Reasonable agree- 
ments of the form factor- and collinear limit of Eq. d38l ) with 
the data fits are found for the parameters m q = 0.834 GeV, 
m s = 0.632 GeV, A = 0.067 GeV, A = 0.448 GeV, n = 0.971, 
g n = 3.604. Details of the fitting procedure for this and an anal- 
ogous calculation for the Sivers function will be presented in a 
future publication 1 80] . 

With the predicted GPD H* and the lensing function 

I'(x,Bt) = b l T j\bj\I(x, \bj\) as input we use Q to give a pre- 
diction for the valence contribution to the first kj -moment of 



the pion Boer-Mulders function, 

db T b 2 T I(x,b T )—<H*(x,b 2 T ). (39) 
o obj 

Numerical results for m 2 h^ l) (x) are presented in Fig. [3] for a 
U(l), SU(2) and SU(3) gauge theory. One observes that all re- 
sults are negative which reflects the sign of the lensing function. 
It was argued in Ref. B66I1 that a negative sign of the lensing 
functions indicates attractive FSIs. We find that this is valid in 
an Abelian perturbative model as well as our non-perturbative 
model in a non-Abelian gauge theory. The magnitude of the 
SU(3) result is about 0.01, while the SU(2) result and U(l) 
result are smaller. One observes a growth of the pion Boer- 
Mulders function with N c . Similar growth was also predicted 
by a model-independent large N c analysis for the nucleon [92], 
though with different leading order behavior. A phenomeno- 
logical calculation of the ratio kxhf w (x,kx)/m„F*(x,k x ) was 
carried out in a perturbative, Abelian, one gluon exchange ap- 
proximation and used to estimate the cos 2<p AA in n~p Drell- 
Yan scattering [14 111 . Similar effects were seen as compared to 
pp Drell Yan scattering Jgl H IH It will be useful to 
study the dependence of the cos 20 AA on color degrees of 
freedom from this non-perturbative approach. So far the pion 
Boer-Mulders function is an unknown function but may be ac- 
cessible from a proposed pion-proton Drell- Yan experiment by 
the COMPASS collaboration. If a pion Boer-Mulders function 
is extracted from such an experiment our analysis can be used 
to quantitatively test the GPD - TMD relation (0. As a compar- 
ison, an extraction of another T-odd parton distribution, the pro- 
ton Sivers function f^J; \ from SIDIS data measured at HER- 
MES and COMPASS J^B reveals an effect of the magnitude 
of about 0.04. A similar calculation using eikonal methods for 
the proton Sivers function will be reported elsewhere [80]. 



7. Conclusions 

In this paper we examined the FSIs of an active quark in a 
pion which are essential to generate a non-vanishing chirally- 
odd and (naive) T-odd parton distribution i.e. the Boer-Mulders 
function. We considered a pion in a valence quark configu- 
ration and worked in a spectator framework. The FSIs then 
were modeled by a non-perturbative 2 — > 2 scattering ampli- 
tude which we calculated using eikonal methods. This treat- 
ment sums up all soft gluons between re-scattered quark and 
antiquark while taking into account color degrees of freedom 
and leads to a more complete description of FSIs as compared 
to calculations in the perturbative Abelian one gluon exchange 
approximation. We find that under the kinematical conditions 
of soft gluon exchange for FSIs, the Boer-Mulders function can 
be split into FSIs modeled by eikonal methods and a spatial 
distribution of quarks in a plane transverse to the direction of 
motion. This spatial distribution is described by a chirally-odd 
pion impact parameter GPD which we calculate in the specta- 
tor model. Together, both effects, i.e. FSI and spatial distortion 
give a prediction for the first moment of the Boer-Mulders func- 
tion that can be tested in pion-proton Drell- Yan experiments. 
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